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Nuclear electromagnetic currents derived in a chiral-effective-field-theory framework in- 
cluding explicit nucleons, A isobars, and pions up to N-^LO, i.e. ignoring loop correc- 
tions, are used in a study of neutron radiative captures on protons and deuterons at 
thermal energies, and of A=2 and 3 nuclei magnetic moments. With the strengths of 
the A-excitation currents determined to reproduce the n-p cross section and isovector 
combination of the trinucleon magnetic moments, we find that the cross section and 
photon circular polarization parameter, measured respectively in n-d and ft-d processes, 
are significantly underpredicted by theory. 

PACS Nos.: 12.39.Fe, 13.40.-f, 25.10.+S, 25.40.Lw 



1. Introduction 

Nuclear electromagnetic currents have recently been derived in a chiral-effective- 
field-theory (xEFT) framework including nucleons, A isobars, and pions Formal 
expressions up to one loop have been obtained in time-ordered perturbation theory 
with non-rel ativist ic Hamiltonians constructed from the chiral Lagrangian formula- 
tion of Refs. I ^ l *^!^! Thus, the study in Ref. ^ is similar to the work of Park et al. 
albeit it uses a different formalism. 

The present talk is a much abridged summary of Ref. The currents up to 
next-to-next-to-leading order (N^LO), that is ignoring loop corrections which enter 
at N'^LO, are used to calculate the magnetic moments of A=2 and 3 nuclei, and the 
thermal neutron radiative captures on protons and deuterons. Realistic two- and 
three-nucleon (for A=3) potentials are used to generate the bound and continuum 
wave functions. To have an estimate of the model dependence arising from short- 
range phenomena, the variation of the predictions is studied as function of the cutoff 
parameter, which is used to regularize the two-body operators, as well as function 
of the input potentials — either the Argonne uis (AV18) or CD-Bonn (CDB) in 
combination with respectively the Urbana IX El or Urbana IX* ^ — used to generate 
the wave functions (the AV18 and CDB have rather different short-range behaviors). 

We find that the N^LO calculations do not provide a satisfactory description 
of the experimental data, particularly for the ^II(n, 7)^11 process. It remains an 
interesting question whether N^LO corrections will resolve the present discrepancies 
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Fig. 1. Diagrams illustrating LO [panel a)], NLO [panels b)-c)], and N^LO [panels d)-f)] currents. 
Nucleons, A isobars, pions, and photons are denoted by solid, thick-solid, dashed, and wavy lines, 
respectively. 

between theory and experiment. 



2. Currents up to N^LO 

The currents up to N^LO are illustrated by the diagrams in Fig. [U where we show 
only one of the possible time orderings. The LO and NLO currents, panels a) and 
b)-c), are well known and will not be given here. At N^LO there is a contribution 
originating from (Q/M)'^ corrections to the LO (one-body) current {Q is the low 
momentum scale and M ~ 1 GeV is the typical hadronic mass scale); it reads ^: 
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where the momenta and are defined as = p'^ — p,; and K,; = (p'^ + Pi)/2, 
q is the photon momentum, and 



Cat,* = (1 + Tj,^)/2 , KN,i ^ (us + nvn^z)/2 , fJ-NA = eN,2 + fJ.N,i , (2) 

with fi;5=-0.12 n.m and Ky=3.706 n.m.. 

The N^LO contributions, represented by diagrams e)-f) and involving A-isobar 
excitation, read in the static limit 21 

e Ca <T2 • k2 
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where gA and are the nucleon axial coupling and pion decay constants, Ca = 
M* 9a ^a/F^ and Ca^ = /i* Dt, h-A and /i* are the iV-A transition axial coupling 
constant and magnetic moment, and A denotes the mass difference to a — m^r. 

The configuration-space versions of the NLO and N^LO operators have l/r^ and 
1 /r^ singularities (r is the interparticle separation) , which need to be regularized in 
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order to avoid divergencies in the matrix elements of these operators between nuclear 
wave functions. We adopt a simple regularization procedure, i.e. a momentum-space 
cutoff. While its precise functional form is arbitrary, the choice made here of a 
Gaussian cutoff fmiction, Ca(p) — e~'^P/^^ , with the parameter A < M, is merely 
dictated by convenience, since it leads to analytical expressions for the Fourier 
transforms It is expected that this arbitrariness be of little relevance, since the 
dependence of theoretical predictions on variations in the cutoff is (or should be, see 
next section) largely removed by a renormalization of the theory free parameters, 
which are fixed by reproducing a given set of observables. 

3. Results and Conclusions 

At N'^LO, the only isoscalar terms are from the (one-body) LO and N^LO-RC oper- 
ators, which are independent of the cutoff A. In Tables[l]we list their contributions 

Table 1. Contributions in n.m. to the deuteron (columns 
labeled AV18 and CDB) and '^He/^H isoscalar combination 
(columns labeled AV18/UIX and CDB/UIX*) magnetic mo- 
ments. The experimental values are 0.8574 n.m. and 0.4257 
n.m., respectively. 



AV18 


CDB 


AV18/UIX 


CDB/UIX* 


LO 0.8469 


0.8521 


0.4104 


0.4183 


N^LO-RC -0.0082 


-0.0080 


-0.0045 


-0.0052 



to the deuteron magnetic moment and isoscalar combination of the ■^Hc and •^H 
magnetic moments. The N^LO-RC correction is (in magnitude) about 1% of the 
LO contribution but of opposite sign, so that its inclusion increases the difference 
between the measured and calculated values. As a result the experimental deuteron 
and trinucleon isoscalar magnetic moments arc underpredicted by theory at the 
(1.6-2.1)% and (3.0-4.7)% levels, respectively, depending on whether the CDB and 
CDB/UIX* or AVIS and AV18/UIX combinations are adopted in the A^2 and A=3 
calculations. We note that a recent calculation of these same observables ^1^, based 
on variational Monte Carlo (VMC) wave functions corresponding to the AV18/UIX 
Hamiltonian, finds the magnitude of the N^LO-RC correction somewhat smaller in 
A=2 (-0.0069 n.m.) and significantly larger in A—i (-0.012 n.m.) than obtained 
here. However, the expression for the magnetic dipole operator is different from 
that resulting from Eq. Ilj, and the VMC wave functions are less accurate than the 
hyperspherical harmonics (HH) wave functions used in this work. 

In the isovector sector, the NLO current involves the combination qaIFt^., for 
which we adopt the value (m^r 5A/^7r)^/(47r)=0.075 as inferred from an analysis of 
nucleon-nucleon elastic scattering datal^^J. In the N^LO currents, the parameters 
Ca and Ca^ are determined by reproducing the n-p radiative capture cross section 
and ■^He/'^H isovector magnetic moment, respectively. We note that the Ac current 
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in Eq. (j4]) gives no contribution in the n-p capture. It contributes in three-body 
matrix elements only because in the configuration-space version of this operator 
the (5-function is replaced by a finite width Gaussian It is for this reason that 
one can interpret the contributions resulting from the A^, current as representing 
corrections beyond N'^LO. 

Results for the isovector combination fiy of the trinucleon magnetic moments 
(without inclusion of the Ac current contribution) are presented in Table [2l The 
NLO contribution calculated in Ref. ^ with VMC wave functions and a cutoff of 

Table 2. Contributions in units of n.m. to the isovector combination 
of the trinucleon magnetic moments, obtained with the AV18/UIX and 
CDB/UIX* Hamiltonian models and cutoff values in the range 500—800 
MeV. The LO and N^LO-RC contributions are cutoff independent. The 
experimental value is -2.553 n.m.. 







AV18/UIX 






CDB/UIX* 




A (McV) 


500 


600 


800 


500 


600 


800 


LO 


-2.159 


-2.159 


-2.159 


-2.180 


-2.180 


-2.180 


NLO 


-0.156 


-0.197 


-0.238 


-0.113 


-0.156 


-0.200 


N^LO-RC 


-(-0.029 


+0.029 


+0.029 


+0.024 


+0.024 


+0.024 


N^LO-A 


-0.258 


-0.253 


-0.250 


-0.205 


-0.202 


-0.200 


Sum 


-2.544 


-2.580 


-2.618 


-2.474 


-2.514 


-2.556 



600 MeV is -0.205 n.m., which is 4% larger than obtained here. Of course, the 
parameter Ca^ is adjusted to reproduce, as function of A, the fiy experimental 
value. 

Predictions for the cross section ax and photon circular polarization parameter 
Rc measured in the reaction ^II(n, 7)^^11 (with unpolarized and polarized neutrons, 
respectively) are presented in Table [3l At thermal energies this process proceeds 
through S-wave capture predominantly via magnetic dipole transitions from the 
initial doublet J=l/2 and quartet J=3/2 n-d scattering states. In addition, there is 
a sm all c ontribution due to an electric quadrupole transition from the initial quartet 
statellSl. At N^LO the cross section is underpredicted by theory by (11-38)% as the 

Table 3. Cumulative contributions to the cross section cry (in mb) and 
photon polarization parameter Rc of the reaction ^H(n,7)^H at thermal 
energies, obtained with the AV18/UIX Hamiltonian model and cutoff values 
in the range 500-800 MeV. The experi men tal values for (jy and Rc are 
respectively (0.508 + 0.015) mb from Ref. El and -0.42 + 0.03. from Ref.HSl. 













Rc 




A (MeV) 


500 


600 


800 


500 


600 


800 


LO 


0.229 


0.229 


0.229 


-0.060 


-0.060 


-0.060 


LO+NLO 


0.272 


0.260 


0.243 


-0.218 


-0.182 


-0.123 


LO+ • ■ ■ +N2LO-RC 


0.252 


0.241 


0.226 


-0.152 


-0.109 


-0.041 


LOH hN^LO-A 


0.438 


0.416 


0.389 


-0.432 


-0.418 


-0.397 


LOH hN^LO-Ac 


0.450 


0.382 


0.315 


-0.437 


-0.398 


-0.331 
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cutoff is increased from 500 MeV to 800 MeV. This rather drastic cutoff dependence 
is mostly due to the contribution of the N^LO-Aj, current. Indeed removing it leads 
to a much weaker variation of the cross section — roughly ±5% about the value 
obtained with A = 600 MeV (next to last row of Table [3|) . It will be interesting 
to see to what extent, if any, loop corrections at N'^LO will improve the present 
predictions, and in particular reduce the cutoff dependence. 

The photon polarization parameter is very sensitive to contributions of NLO and 
N^LO currents, which produce more than a sixfold increase, in absolute value, of 
the LO result, and bring it into much closer agreement with the measured value. All 
results listed in Table[3]for Rc (and ctt) include the small 644 RME, although it only 
has a significant effect for the LO prediction (_Rc=^0.060 versus -0.072 depending 
on whether 644 is retained or not). 

We conclude by summarizing our results. Up to N'^LO, the only isoscalar terms 
are those generated in a non-relativistic expansion of the one-body current, and 
provide a (cutoff- independent) 1% correction — relative to LO — to the deuteron and 
isoscalar combination of the trinucleon magnetic moments. This correction is of 
opposite sign to the LO contribution, and therefore increases the underprediction 
of the corresponding experimental values from (0.9 ± 0.3)% for the deuteron and 
(2.7 ±0.9)% for the trinucleons at LO to, respectively, (1.9 ±0.3)% and (3.8 ±0.8)% 
at N^LO. The spread reflects differences in the short-range behavior of the AV18 and 
CDB potentials, in particular the weaker tensor components of the latter relative 
to the former in this range. 

At NLO, isovector terms arise from the pion seagull and in-flight contributions, 
while at N^LO, in addition to the relativistic corrections mentioned above, isovector 
terms due to A-isobar excitation are also obtained. The parameters Ca and Ca^ 
of the N^LO two-body A-excitation currents have been determined, as functions of 
the cutoff A and for the Hamiltonian model of interest, by reproducing the cross 
section for the n-p radiative capture at thermal neutron energies and the isovector 
combination of the trinucleon magnetic moments. This current has then been used 
to make predictions — with the AV18/UIX model only, since HH continuum wave 
functions are not yet available for the CDB/UIX* model — for the cross section 
fjT and photon circular polarization parameter Rc measured in the capture of, 
respectively, unpolarized and polarized neutrons by deuterons. The experimental 
gt {\Rc\) is found to be underestimated by 11% (overestimated by 4%) for A=500 
MeV and 38% (underestimated by 21%) for A=800 MeV. 

The results display a significant cutoff dependence, particularly so for the N^LO 
contributions associated with A isobar degrees of freedom. Indeed these contribu- 
tions are much larger than those at NLO. This is partly due to the fact that the 
two NLO (pion seagull and in-flight) terms interfere destructively. For example, 
the seagull (in-flight) contributions to doublet 77122 and quartet TO44 Ml matrix 
elements, in units of ivc?^^ and for A=500 MeV, are respectively -9.1 (+6.5) and 
-0.8 (+0.6). As a result cjt = 0.425 mb and Rc = -0.425 at LO+NLO (seaguU 
only), which should be compared to (Jt = 0.272 mb and Rc = —0.218 at LO+NLO 
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(seagull+in-flight) from the second row of Table[3l The relatively large A-excitation 
contributions also point to the need for including loop corrections at N^LO, which 
these N^LO currents, because of the procedure adopted here to determine their 
strength, are implicitly making up for. 

The next stage in the research program we have undertaken is to incorporate the 
N'^LO operators derived in Ref. ^into the calculations of the captures and magnetic 
moments involving light nuclei (with mass number A < 8), and indeed to extend 
these calculations to also include p-d capture at energies up to a few MeV's, and 
possibly four-nucleon processes, in particular '^He(n, 7)^He at thermal energies. Of 
course, at N'^LO three-body currents also occur, and will need to be derived. Work 
along these lines is being pursued vigorously. 
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